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Abstract 

The metric of a Schwarzschild solution in brane induced gravity in five 
dimensions is studied. We find a nonperturbative solution for which an exact 
expression on the brane is obtained. We also find a linearized solution in the 
bulk and argue that a nonsingular exact solution in the entire space should 
exist. The exact solution on the brane is highly nontrivial as it interpolates 
between different distance scales. This part of the metric is enough to deduce 
an important property - the ADM mass of the solution is suppressed com- 
pared to the bare mass of a static source. This screening of the mass is due 
to nonlinear interactions which give rise to a nonzero curvature outside the 
source. The curvature extends away from the source to a certain macroscopic 
distance that coincides with the would-be strong interaction scale. The very 
same curvature shields the source from strong coupling effects. The four di- 
mensional law of gravity, including the correct tensorial structure, is recovered 
at observable distances. We find that the solution has no vDVZ discontinuity 
and show that the gravitational field on the brane is always weak, in spite of 
the fact that the solution is nonperturbative. 



1 Introduction, Discussions, and Summary 



Exact static solutions in models of gravity carry a great deal of information on 
the gravitational theories themselves. Hence, finding these solutions in models that 
modify gravity at large distances is an important and interesting task. In the present 
work we will study the Schwarzschild solution in the DGP model U where gravi- 
tational interactions are modified at large cosmological distances. It is complicated 
to find this solution since even at scales much larger than the Schwarzschild radius 
of a source, full nonlinear treatement is required [2]. The first approximate solution 
was obtained in Ref . |2| and subsequently by the authors of Refs. jU IHl CI ■ The 
solution should interpolate between very different distance scales. These scales are: 
the 4D gravitational radius of the source of a mass M, tm = 2G]sfM, the large 
distance crossover scale Vc ~ 10^*^ cm, and an intermediate scale, first discovered by 
Vainshtein in massive gravity jH], which in the DGP model reads as follows 



This is a scale at which nonlinear interactions in a naive perturbative expansion in 
Gn become comparable with the linear terms (we will discuss below the physical 
meaning of this scale in detail). For a source such as the Sun, the hierarchy of the 
scales is as follows: 



In most of the work, unless stated otherwise, we will be discussing sources that are 
smaller than r^. In Refs. [H 12 El IH ISl IHl 171 the approximate solutions for such 
sources were found in different regions of (0). The main properties of the solution 
can be summarized as follows: 

(a) At distances r ^ rc the 5D Schwarzschield solution with the 5D ADM mass 
M is recovered (throughout this work r stands for a 4D radius). 

(b) For <^ r <^ Tc the potential scales as in the 4D Schwarzschild solution. 
However, relativistic gravity is a tensor-scalar theory that contains the gravitation- 
ally coupled scalar mode (i.e. the tensorial structure is that of a 5D gravitational 
theory which contains extra polarizations). 

(c) For r <^ the theory reproduces the Schwarzschild solution of 4D General 
Relativity (GR) with a good accuracy. 

Perhaps the most important property of the (a-c) solution outlined above is 
the dynamical "self-shielding" mechanism by which the solution protects itself from 
the would-be strong coupling regime 0. Very briefly, the self-shielding can be 
described as follows: the expansion in On breaks down at the scale r ~ r* making 
the perturbative calculations unreliable below this scale. However, exact nonlinear 
solutions of equations of motion - which effectively re-sum the series of classical 
nonlinear graphs - are perfectly sensible well below the scale r^,. Hence, the correct 
way of doing the perturbative calculations is first to find a classical background 
solution of equations of motion and then expand around it. 




(1) 
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In the present work we attempt to find exactly the Schwarzschild solution in 
the DGP model. We managed to find explicitly only a 4D part of the metric. 
This exact result, combined with reasonable boundary conditions in the bulk, is 
sufficient to determine unambiguously a number of crucial properties of the solution. 
First we confirm the existence of the scale r^, - this scale enters manifestly our 
exact solution. We also confirm that the self-shielding mechanism outlined above 
takes place. Furthermore, we obtain deeper insight into the dynamics of the self- 
shielding, which, to the best of our knowledge, has not been emphasized so far in 
the literature: the self-shielding effect takes place because a source creates a nonzero 
scalar curvature that extends outside the source to a distance r ~ r^. This curvature 
suppresses nonlinear interactions that otherwise would become strong at the scale 
below r^,. On the other hand, we also find that some of the physical properties of 
our solution differ from those in (a-c). Our solution has the following main features: 

(A) For r ^ Tc, like in (a), one recovers the 5D Schwarzschild solution, however 
unlike in (a), the new solution has the screened 5D ADM mass 



The screened mass is suppressed compared to the bare mass M. Therefore, the new 
solution is energetically favorable over the (a-c) solution. 

(B) For r=K ^ r ^ Tc one can think of the solution as being a four- dimensional 
one with an r-dependent decreasing mass M(r) ~ t^.tm/t. Alternatively, one can 
simply think of the solution just approaching very fast the 5D Schwarzschild metric 
with the screened mass (jH}, i.e., approaching the asymptotic behavior of (A). 

(C) For r <^ r* the results agree with those of (c) with a good accuracy. 

The (a-c) and (A-C) solutions both asymptote to Minkowski space at infinity. 
However, the way they approach the fiat space is different because of the difference 
in their 5D ADM masses. The (A-C) solution, or any of its parts, cannot be obtained 
in linearized theory, it is a nonperturbative solution at any distance scale. Since the 
mass of the (a-c) solution is larger than the mass of the (A-C) solution, we would 
expect that the heavier solution will eventually decay into the light one, unless there 
are some topological arguments preventing this decay. 

The above findings suggest that Minkowski space, although globally stable in 
the DGP model, is locally unstable in the following sense. A static source placed on 
an empty brane creates a nonzero scalar curvature around it. For a source of the 
size ^ r=K this curvature extends to a distance ~ r*. Above this scale the solution 
asymptotes very quickly to 5D Minkowski space. More intuitively, a static source 
distorts a brane medium around it creating a potential well, and the distortion 
extends to a distance r ~ r^,. Since is much bigger than the size of the source 
itself, we can interpret this phenomenon as a local-instability of the fiat space. This 
local-instability, however, has not been seen in the linearized theory It should 
emerge, therefore, in nonlinear interactions and should disappear when the scale 




(3) 
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is taken to zero^. 

It is remarkable that the distance scale to which the local-instability extends, 
coincides with the scale r^, at which the naive perturbative expansion in G^v breaks 
down. Therefore, by creating a scalar curvature that extends to r ~ r*, the source 
shields itself from a would-be strong coupling regime that could otherwise appear at 
distances r ^ [2]: (i) The coupling of a phenomenologically dangerous extra scalar 
polarization of a 5D graviton to 4D matter gets suppressed at distances r ^ r^, due to 
the curvature effects. This is similar to the suppression of the extra polarization of a 
massive graviton on the AdS background |TT1^]. Indeed, in our case the curvature 
created by the source, although coordinate dependent, has the definite sign that 
coincides with the sign of the AdS curvature. As a result, the model approximates 
with a high accuracy the Einstein gravity at r -C with potentially observable 
small deviations fSl Ej (see comments below), (ii) The self-coupling of the extra 
polarizations of a graviton, which on a fiat background leads to the breakdown of 
a perturbative expansion and to the strong coupling problem, gets now suppressed 
at distances ^ r* by the scalar curvature created by the source. This is also similar 
to the suppression of the self-coupling of the massive graviton polarizations on the 
AdS background [MITH] . 

The main properties of the classical solution described above seem to be universal 
and should be expected to hold in models that modify gravity at large distances. To 
see the viability of this argument let us look at the 4D part of the Einstein equations 
of the DGP model 

Gf,u + [K^y - g^yK] = Svr Gn T^^ , (4) 

where G^^ is a four- dimensional Einstein tensor, T^y is the matter stress-tensor, 
K^u{g) is a symmetric tensor of the brane extrinsic curvature and K denotes its 
trace (in other models, e.g., in massive gravity, the extrinsic curvature part will 
be replaced by the mass term). Consider a localized source, a star for instance. 
Outside the source the r.h.s. of (jU is zero. However, the term [K^^ — g^uK] need 
not be zero. This would lead to a nonzero Ricci tensor and scalar. This is unlike 
the Einstein gravity where only the Riemann tensor components are nonvanishing. 
The curvature that is produced away from the source, however, is small since it is 
proportional to the strength of the source itself multiplied by additional suppression 
factors proportional to powers of rric. According to our exact solution the nonzero 
curvature extends effectively to distances of the order r^:, but it is sub-dominant to 
the standard 4D Schwarzschild contribution to the Kretschman scalar, except in the 
region around r ~ r^,, where the two curvature invariants are roughly of the same 
order ~ ml = ~ (10~^^ GeV)^ (see Fig. 2 and detailed discussions in section 
3). One important property of the solution is that the gravitational field is weak 
everywhere outside the source. Nevertheless, the solution is nonperturbative and an 

^The latter assertion is valid since the (A-C) solution, as we will see, is regular in the rric 
limit where it turns to the conventional 4D solution, i.e., it has no vDVZ discontinuity pill(J|. 
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expansion in Gn cannot be used to recover the solution even at very large distances 

In the present paper we are primarily concerned with classical sources. Nev- 
ertheless, we would like to comment as well on dynamics of "quantum" sources, 
such as gravitons. Consider the following academic set-up: a toy world in which 
there is no matter, radiation and/or any classical sources of gravity - only gravitons 
propagate and interact with each other in this world. Because of the very same tri- 
linear vertex diagram that leads to the breakdown of the expansion for classical 
sources (see Ref. j2]), the self-interactions of the gravitons will become important 
at lower energy scale than they would in the Einstein theory. The corresponding 
breakdown scale is the scale (0) adopted to a quantum source with tm = 1/Mpi, 
that is ~ (r^/Mpi)^/^ (see also Ref. [12] that obtains a somewhat different 
scale). In this set-up the graviton loop diagrams could in principle generate higher 
derivative operators that are suppressed by the low scale. A theory with such high- 
derivative operators would not be predictive at distances below A^^ ~ 1000 km or 
so. 

However, there are two sets of arguments suggesting that the above difficulty 
might well be unimportant for the description of a real world which, on top of the 
gravitons, is inhabited by planets, stars, galaxies etc. We start with the arguments 
of Ref. j7j . This work takes a point of view that A^ is a true ultraviolet (UV) cutoff 
of the theory in a sense that at this scale some new quantum gravity degrees of 
freedom should be introduced in the model. Nevertheless, as was discussed in detail 
in Ref. this should not be dangerous if one considers a realistic setup in which 
mater is introduced into the theory. For instance, consider the effect of introducing 
the classical gravitatonal field of the Sun. Because of the gravitational background 
of the Sun, the UV cutoff of the theory becomes a coordinate dependent quantity 
Aq(x). This cutoff grows closer to the source where its gravitational field becomes 
more and more pronounced, hence, increasing the value of the effective UV cutoff. 
In this approach the authors of Ref. [7j managed to find a minimal required set 
of higher dimensional operators that are closed with respect to the renormalization 
group fiow. Because of the re-summation of the large classical nonlinear effects these 
operators are effectively suppressed by the coordinate dependent scale Aq(x). If so, 
the new UV physics won't manifest itself in any measurements j7]. 

Putting all this on a bit more general ground, one should define the model in 
an external background field. That is, in the action and the partition function of 
the model the metric splits into two parts Qf^u = fi'^i + g'jiu^ where g'^^ stands for 
the classical background metric and g'j^^, denotes the quantum fiuctuations about 
that metric. The classical part satisfies the classical equations of motion with given 
classical gravitational sources such as matter, radiation, planets, stars, galaxies, 
etc... Then, the effective UV cutoff for quantum fluctuations at any given point in 
space-time is a function of the background metric. For a realistic setup this effective 
cutoff is high enough to render the model consistent with observations. 

We find the above logic useful and viable. We also think that the algorithm of 
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Ref. j7| might be the most convenient one for practical calculations. Nevertheless, 
there could exist deeper dynamical phenomena to the discussions of which we turn 
right now. Although our arguments below parallel in a certain respect those of 
Ref. [7], there is a conceptual difference on the main issue. Our view, that we will 
try to substantiate in subsequent works, is that the scale Ag is not a UV scale of 
the model in the sense that some new quantum gravity degrees of freedom should 
be entering at that scale. We think that all what's needed to go above the scale 
Ag is already in the model, and that this is just a matter of technical difficulty 
of non-perturbative calculations (or, in other words, is a matter of difficulty of 
summing up loop diagrams). The re-summation could in principle cure problems at 
the loop level as well. At this end, we do not see a reason why the self- shielding 
mechanism outlined above should not be operative for "quantum" sources too. The 
very same local-instability of Minkowski space should manifest itself in nonlinear 
interactions of quantum sources, e.g., gravitons. The local-instability scale in that 
case is Ag. Hence, we would expect that a quantum source creates a curvature 
around it that extends to the distances of the order of ~ (r^ 

/Mpi)i/3 ^ 1000 km, 

and doing so it self-shields itself from the strong coupling regime. If this is so, 
then the problem of loop calculations boils down to the problem of defining correct 
variables w.r.t. which the perturbative expansion should be performed. In this case 
the field decomposition should take the form: g^i, = + (7^^, where stands for 
a nonperturbative background metric created by a "quantum" source. Similar in 
spirit arguments using a toy model were given by Dvali in Ref. j l7J . 

In this work we solve exactly for the 4D part of the Einstein equations of the 
DGP model. Furthermore, we study the bulk metric as far as we can. Here we 
do not have an exact solution. Nevertheless, a number of important and reliable 
properties can be deduced. Using the analytic continuation and taking advantage of 
the fact that the bulk asymptotes to Minkowski space, we argue that a nonsingular 
solution that matches our brane solution should exist in the bulk. Furthermore, we 
find large- distance asymptotes of the bulk solution. It is important to point out 
that irrespective of the form of the bulk solution (as long as it is nonsingular), we 
are able to deduce the properties (A-C) of the system. 

Finally, we would like to make two important comments. First, the DGP model 
possesses two branches of solutions that are distinguished from each other by the 
bulk boundary conditions. These two branches are disconnected. In this work we 
concentrate primarily on the Schwarzschild solution of the so-called conventional 
branch on which the brane and the bulk asymptote to Minkowski space at infin- 
ity. However, the second, the so-called "self-accelerated" branch [18j is extremely 
interesting as it can be used to describe the accelerated expansion of the Universe 
without introducing dark energy ^^l- In the present work we also find an exact 
brane metric for a Schwarzschild source on the self-accelerated branch. However, 
because the asymptotic behavior of the solution on this branch is not Minkowski we 
are not able to argue for the existence of a nonsingular bulk solution. On the other 
hand, we do not see any physical reason why this solution should not exist in the 
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bulk as well. This branch will be discussed in detail elsewhere. Second, it is inter- 
esting to note that the linearized analysis of the DGP model in dimensions six and 
higher [201 , as well as certain modifications of the five-dimensional model 123 
show no sign of breakdown of perturbation theory and strong nonlinear effects. It 
is left for future work to understand more deeply the interconnections between all 
these approaches. 

The paper is organized as follow. In Section 2 we set the action and equations of 
motion of the DGP model. In section 3 we give a qualitative description of the main 
new properties of our solution. In section 4 we give exact solutions for the metric 
and extrinsic curvature on the brane. In section 5 we discuss the absence of the 
vDVZ discontinuity and in section 6 we comment on distinctions between the Gjv 
and the weak-field expansions in the model. In section 7 we use the ADM formalism 
to argue that the solution for the metric and extrinsic curvature can be smoothly 
continued into the bulk space. Concluding remarks are given in Sect. 8. 



2 The setup 

We consider the action of the DGP model [11: 

S = Ml j d^xy/^R + Ml j d^x^f^R . (5) 

Here, the (4 + 1) coordinates are x^^ = {x^, y), = 0, . . . , 3 and g and R are the 
determinant and curvature of the 5 dimensional metric Qmn, while g and R are the 
determinant and curvatures of the 4 dimensional metric g^u = g^uix^, y = 0). 
The (fiu) and (yM) equations of motion are respectively, 

y-g G^u{x)6{y) + G^i,{x,y) = , (6) 

GyAiix^y) = 0, (7) 

where rric = 2M^/Mp is the inverse of the crossover scale (the Gibbons- Hawking 
surface term [23 that guarantees the correct Einstein equations ()6|7|) is implied in 
©)• 

We will study the analog of the Schwarzschild solution in this setup. Thus, we 
consider the most general static metric with spherical symmetry on the brane and 
with Z2 symmetric line element: 

ds'^ = -e^dt^ + e^dr^ + ef'r^dn^ + 7 drdy + e^dy^ , (8) 

where i>, A, /i, 7, a are functions of r = ^J^xFg^il^ and y. The Z2 symmetry across 
the brane implies that 7 is an odd function of y while the rest are even. 

The jump conditions on the y derivatives of the warp factors across the brane 
(at y = 0) implied by ^ give^ 



2e-^+t + 2e"^ 



1 3 
-1 + rA,. - 3r/i,. + - -^l^r^l^r - ^V' 



In the second line of (|5J, (|10() and Hll() the functions are evaluated at y = 0+. 
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corresponding to the (tt), (rr) and (99) components. 

We have not made use of any gauge freedom up to this point. A convenient 
choice is to set /i(r, = by rescahng r — >• r exp(— /i/2). Moreover, we set z/(r, y) = 
—X{r,y) by transformation of the y coordinate. The resuhing hne element is (note 
that this gauge can be reached on the brane and in the bulk only because we allowed 
an off-diagonal term in the metric) 



where A, 7, a are functions of r and y. Our choice is such that the brane is 
not bent in this coordinate system and is located at y = 0. A more conventional 
diagonal coordinate system can be obtained by a coordinate redefinition after which 
the interval reads 



Here the functions u and f3 are related to A, 7, and a. In the z, r coordinate system 
the brane is bent. Typically in brane- world models the 4D part of the Einstein 
equations are not closed. Hence, the induced metric on a brane cannot be deter- 
mined without some input from the bulk equations, and/or without making certain 
assumptions about the induced metric itself. This would also be true in our case. 
However, in the gauge ()12|) . we find a subset of the Einstein equations that can be 
closed for the function A. As a result, A can be found exactly on the brane. Al- 
though the knowledge of A alone is not enough to describe the whole gravitational 
dynamics on the brane ( for instance, this is not enough for the description of 5D 
matter geodesies at short distances since transverse derivatives of the metric are 
also entering the 5D geodesic equations) nevertheless, combining the knowledge of 
A with the asymptotic behavior of the other functions in (fT^ that we can also ob- 
tain unambiguously, is enough to deduce the properties (A-C) of the Introduction. 
Hence, these properties are "exact". 

Once a source is placed on the worldvolume, the brane produces a nonzero extrin- 
sic curvature. As a result of this a nonzero 4D intrinsic curvature is also produced. 
We will discuss these issues in detail in the next section. 



ds^ = -e'^dt^ + e^dr^ + r'^d^f + 7 drdy + e'^rfy^ . 



(12) 



ds 



(13) 
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From (|7j) we also derive the following equation (in the gauge 

+ (e^+'^ - 7^) (472 + (4r7 - r^e^Xy - r^^Xr)Xy + 2r^^Xry) 

+ e'^(27^ - -far - lXr){Ar-f + r^e^Xy - r^-fXr) . (14) 

3 Structure of the solution 

In the next section we solve the above system of equations exactly on the brane and 
obtain a perturbative solution in the bulk. The fact that the exact solution could 
be found on the brane is nontrivial. However, the solution can only be written in 
an implicit form from which the extraction of useful information requires additional 
efforts. To simplify the reading of the paper we summarize certain nontrivial prop- 
erties of the solution in the present section. The solution itself will be derived in 
the next section. 



3.1 Solution on the brane 

In this subsection we discuss the properties of the solutions on the brane, i.e., at 
y = 0. We find certain similarities, as well as differences, in the 4D part of our 
solution with the anti-de Sitter-Schwarzschild (AdSS) solution of conventional 4D 
General Relativity (GR) with a small negative cosmological constant A 

Consider 4D GR with the cosmological constant A = — 3m^. Furthermore, con- 
sider a static source of mass M (a star) and a Schwarzschild radius tm = 2GnM in 
this space. In the static coordinate system the AdSS solution takes the form 

ds' = -(1 - + ml A de + ^ + rUnl. (15) 



(1-^ 



This coordinate system covers the AdSS solution in the interval tm < r < Vf. = m~^. 
The following properties of the AdSS solution will be contrasted to our solution, 
(i) In the interval tm < r < there is a new distance scale r^, exhibited by (jisp 

^MrlY^^ ■ (16) 



The physical meaning of this scale is as follows. For r < the Newtonian potential 
r^/r in ()15|) dominates over the term m^r^, while for r > r^, the term m^r^ overcomes 
the Newtonian term. Hence, is a scale at which the Newtonian and the m^r^ 
terms are equal. This can also be expressed in terms of invariants. Let us define the 
Kretschman Scalar (KS) 



Rk ^ J{Rits)\ (17) 
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where R^'fl^s is a Riemann tensor of the Schwarzschild part of the solution (i.e., of the 
part that survives in the hmit). We compare the KS with the background 

curvature due to the cosmological constant 

\Ra\ = 12ml. (18) 

We get 

Rk:^\Ra\ for r<r*; Rk <^ \Ra\ for < r < . (19) 

Therefore, r^, is a scale at which Rx — \Ra\- For r^f ^ r ^ r^, the corrections 
due to the background curvature are small and the solution is dominated by the 
Schwarzschild metric, while for r^, -C r -C Tc the background curvature terms are 
bigger that the Schwarzschild terms, both of them still being smaller than 1. 

(ii) At r 3> the Schwarzschild part becomes irrelevant compared to the AdSS 
part. 

We will show below that our solution has some of the properties described in (i), 
however, unlike (ii), it behaves as 5D Schwarzschild solution at large distances. 
The 4D part of our solution (i.e. the solution at y = 0) for r <^ takes the form 



ds 



2| 



\y=0 



•-f + rnyg(r)) + f , . , ^^ + dill (20) 



Like the AdSS solution, the metric (pUj) possesses the r^, scale defined in (fT^ . As 
we will see below this scale has the same physical meaning as in the AdSS case. For 
instance, at r <^ r^, 



1 

,4' 



9ir) - [^j . (21) 

Then, it is straightforward to check that 

^ mlr'^gir) for r <^ r^,; — ^ ~ mlr'^gir) for r ~ . (22) 

Therefore, the corrections become of the order of the tm/t term at around r ~ r^,. 
Moreover, like the AdSS solution, the corrections dominate over r^z/r for r^, <^ 
r <C turning the 4D behavior of the solution into the 5D behavior. The plot 
of the function is given on Fig. 1. As in the AdSS case, the corrections to the 
Schwarzschild solution that are proportional to mc give rise to the four- dimensional 
Ricci curvature Rm^- This is interesting since the curvature is completely due to 
the modification of gravity. However, unlike the AdSS case, this curvature is not a 
constant but depends on r; moreover it also depends on the strength of the source 
itself. The plot of the Ricci curvature is given on Fig. 2. 
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Figure 1: The solid line depicts P{r)/rM (on the vertical axis), wher P is defined in 

as a function of r on the horizontal axis. The dashed line presents the function 0.28r=K/r; 

The value of is set to 1 on this graph. 




Figure 2: The solid line depicts the magnitude of the four-dimensional Ricci scalar cur- 
vature (on the vertical axes) as a function of r on the horizontal axes. The dashed line 
depicts the dependence of the 4D Kretschmann scalar on r. The value of is set to 1 on 
this plot. 
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Similar to the AdSS solution the above properties can be expressed in terms of 
the invariants 

Rr > Rmc for r <^r^ and for r > r^, ; while Rk ~ Rmc for r ^ . (23) 

Unlike in the AdSS solution, however, the curvature Rm^{r) decreases very fast after 
r ^ r*. Hence, the induced curvature Rm^{r) is sub-dominant to Rk everywhere 
except in the neighborhood of the point r ~ r^, where both of these are of the same 
order ~ m^, see Fig. 2. 

Furthermore, unlike the AdSS solution, our solution can be presented in the same 
coordinate system even for r > Tc- This is because there is no horizon at r = rc and 
our solution smoothly turns into the 5D Schwarzschild solution 

ds\^o = - f 1 - J) + + dnl . (24) 

The key property of this solution is that the gravitational radius is rescaled 

tm ^ tm [ — ] > • (25) 

This has an explanation. The gravitational radius grows compared to tm because 
in the 5D regime the gravitational coupling constant grows. However, there is an 
opposite effect as well. In fact, the gravitational radius reduces compared to what 
it would have been in a pure 5D theory with no brane. This is because the effective 
mass of the source Meg, defined as f\j = Mes/M^, gets screened. 

All the above results could be understood as follows. Consider an empty brane 
and an empty bulk. Minkowski space is a solution. Let us localize a static source 
on the brane at r = 0. The Minkowski solution remains globally stable, however, 
the source, no matter how weak, triggers local instability of Minkowski space on 
a brane in the region r < r^,. In this patch Minkowski space is readjusted to a 
curved space. The curvature of the latter depends on the strength of the source, it 
slowly decreases with increasing r but drops fast at r > r^,. For an observer at large 
distances it looks as if the source has polarized the medium (brane) around it. This 
observer measures the screened mass (j25|) which also includes the contributions of 
the curvature. 



3.2 Solution in the bulk 

At large enough distances, i.e. a/PT^ ^ r^, the solution turns into a 5D 
spherically-symmetric Schwarzschild solution: 

^^^'1 rr-^ ~ - f 1 - j'' df + — '^"'\ , + (r^ + y") dnl . (26) 



' M 
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However, the 5D spherical symmetric is only an approximation and does not hold 
for \/r^~+~y^ <^ r^,. In the latter regime the properties of the solution on and off 
the brane are rather different. The pure 5D spherically-symmetric solution 1)2(11] is 
squeezed both on and off the brane but it is more squeezed on the brane then in 
the bulk. Hence, the only symmetry of the solution that is left is the cylindrical 
symmetry. 

4 The solutions 

In this section we solve for A on the brane in the coordinate system (P^. 
From the and (|Tnjl equations we deduce: 

7r = e^K , (27) 

and © and (|TT|) can be rewritten as 



2e^P, = (A,--)7-7r , (28 



A 



/ [K - -n - 7r 

^QX+a _ ^2 V r 

2mcr2 / 1 
/ [K - -n - 7r 



retp,, = ^^£^ ( (A, - -)7 - 7, ) , (29) 



where we have defined 



P(r) = r(l - exp(-A)). (30) 

Using and in the yy equation of motion ((7j) we obtain: 

= Qmy{2Q + rQr) + r'Ql - SrQQr + AQ^ . (31) 

where Q = P^. Let us study this equation. First of all we note that the quadratic 
terms in ()31|) can only be neglected when r is large. This suggests that the naive 
linearized approximation that neglects the quadratic terms is viable only for large 
distances, no matter how weak the source is. Hence, for r ^ we can neglect 
the last three terms on the r.h.s. of (|31|) and then the solution is 

p = !il + Ci , (32) 
r 

where tm and Ci are integration constants. Ci = gives the 5D Schwarzschild 
solution of radius tm- In a similar fashion, for r -C we can neglect the term 
proportional to in ()3T|l . The solutions in this case are 

P = rM + C2r2(2±v^)+i , (33) 
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where tm and C2 are integration constants. In this case, C2 = gives the 4D 
Schwarzschild solution of radius r^. 

We will show that there is an interpolating solution between these two regimes 
(regular branch) together with a second solution that becomes 5D de Sitter Schwarzschild 
at large distances (accelerated branch). In order to find these exact solutions let us 
rewrite as 

= U^- 4(1 + U)U, - 8U{2 + U) , (34) 

where z = ln(r/ro) and U = — 2exp{— 2 z)Q /3mlrl (ro an arbitrary constant). There 
are two solutions that are given implicitly by: 



In 



In 



;i + 3f/ + /) 



f/2(3 + 3f/ + V3/)2v^(-5 - 3f/ + /) 
(-5 -3U + /)(-3 -3U- V3ff^ 



(f/ + 2)2(1 + 3t/ + /) 



8Z + C3 



(35) 
(36) 



where / = VI + 6f/ + 3f/2 and C3 an integration constant. 

Let us first consider the solution For f/ > the left hand side is a decreasing 
function of U . To obtain the large distance (r ^ r^) behavior of the solution we 
consider the limit f/ — > O"*" (2; — > +00) in which (jH^ reads 



In f/ = -Az + C4 + 0{U) , 

that gives ()32|) upon integration. The short distance regime (r <^ r*) 
taking the limit U —>■ +00 {z — >• —00). In this case flH^ is given by 

1 



:i + V3)\n U = -Az + C^ + O (^-) 



(37) 

is found by 



(38) 



that reproduces the asymptotic behavior of ()33|) (with the minus sign). Thus, we 
have found a smooth solution that interpolates between the 4D and 5D Schwarzschild 
solutions on the brane. This corresponds to the regular branch. 

Let us now study the second solution For U < —2 the left hand side of 

fl36|) is an increasing function of U. In the large distance limit U —2^ {z — > +00) 
flHUj) reads 

ln(-f/ - 2) = -Az + C4 + 0{U + 2) , (39) 



that gives upon integration 



(40) 



where we have absorbed the integration constant C4 in the definition of tm- This 
is the 5D de Sitter-Schwarzschild solution of the accelerated branch. The short 
distance behavior is obtained by considering the U —00 {z —00) of (|36|) . This 
gives 



;i + V3)ln(-f/) = -Az + C5 + 0(- 



1 



U 



(41) 
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The above equation gives the following behavior of the potentials for rM ^ r ^ 

z/ = -A ~ - ^ - OMmlr^(^Y^'''\ (42) 
r \ r J 

Below we derive in detail the relation between the mass of the source that deter- 
mines the 4D Schwarzschild radius tm and the "screened" mass that determines the 
5D Schwarzschild behavior (with radius tm) at large distances (r ^ rc). In order 
to do that let us rewrite as 

F{U) = ku , (43) 
by defining the variable u = r^, the constant k = exp(3C3/8)/rQ and 

F(U)= ( ) . (44, 

F{U) is a monotonically decreasing function for f/ > that diverges in the limit 
U —>■ 0+ and vanishes in the limit U +oo. Note that U = — 2P„/m^. On the 
other hand, (|43|) gives U = F~^{ku). Therefore 

- — / Pudu = - / F-\ku)d{ku) , (45) 
JO k Jo 



where the integral on the right hand side is equivalent to the following integral 

F{U)dU ^ 0.43 . (46) 

Imposing the asymptotic behavior of P: P{0) = tm and P(+oo) = on (jlSjl gives 

2krM = OASml (47) 
In turn, the large distance behavior of P: P ~ "rli/r obtained from ()35|) gives 

/ o \ 3/4 

Thus, from ()45p and ()48|) we obtain the exact relation between the 4D and 5D 
Schwarzschild radii: 

TM = 2.6 hi J— . (49) 
V rc 

The 4D Ricci curvature P4 can be readily calculated using the expression: 

R4 = -lml{U, + 4U). (50) 
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It goes to zero in the large r limit (in the regular branch) and grows as one approaches 
the source at small r. Moreover, it gives rise to the properties described in the 
previous section. 

In this gauge we solve for the combination of 7 and a at y = 0^ which takes the 
form 

e~?7 = arctan (—^ — = arctan ^ ^^^° e^(2[/ — f4)^ . (51) 

The latter expression will be used in the next section to check the continuity of the 
solution in the m^, limit. 



5 No vDVZ discontinuity 

The vDVZ discontinuity [HI is an interesting observation that the theory in the 
rric — limit could differ from the one in which rric = is set ab initio (i.e., from 
the Einstein theory). Below we will argue that the vDVZ observation is based on 
arguments that do not hold when the dynamical effects of the mass screening are 
taken into account. 

We will show that our solution is continuous in the rric —>■ limit. We present 
these arguments it two different ways. First let us look at the solution of the previous 
section. In the limit rric (i.e. Vc —>■ 00) we find: 

00 ; Ta/ 00 ; M^s . (52) 

Moreover, the off-diagonal term in the solution behaves as follows: 



7e~"/'|r»r, ^ arctan ( ) < 1 , (53) 

rrirr-^ ' 



which is always small and the regime of the applicability of the above expression 
goes to infinity in the limit rUc 0. Moreover, in the region r -C r^, we get 

_ /2, / rric \ 

76 |r<r. — ^ arctan — g ^0. (54) 

\mc I 

Based on the above findings we conclude that the solution turns into an exactly 4D 
Schwarzschild solution in the rric —>■ limit. The regions where it could deviate from 
the 4D solution, r ^ min{r^,,rc}, go to infinity. 

The vDVZ discontinuity was originally formulated in ^ ^] in terms of a one- 
graviton exchange amplitude. It is instructive to see the loophole in this formalism 
as well. The arguments of jHl CSl go as follows. The momentum space amplitude for 
one graviton exchange between the source of a stress-tensor T^o^^^^e ^ probe T^^ 
is given by 

-j^ ^source rpZ/iu 1 ^source 

Mpi -I- rricP + ie 
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(Here T's denote the trace of the stress-tensors.) The very same amphtude for 
massless 4D gravity is 



-I 'psource rp/^u ]_T^source n^/ 

^' = ^— ^ • ^^6) 

Mpj + le 

Hence, in the hmit mc — Eq. does not reduce to Eq. (jSUj) . This is the vDVZ 
discontinuity. 

Our exact solution suggests that these arguments do not hold in an intricate way. 
Let us start with r ^ r^,. In this region the linearized equations turn into source- 
free equations since the source is only localized at distances r -C r^,. The source-free 
equations can be solved and the solution is a 5D Schwarzschild metric with the 
Schwarzschild radius being a yet unspecified integration constant. This constant 
should either be fixed by matching to the solution at r < r^,, or by calculating the 
ADM mass. In either case the curvature created by the source will also contribute to 
the ADM mass. Therefore, at large distances where the perturbative one-graviton 
exchange is believed to be a good approximation, we have to replace a source T^^"^^ 
with an effective source that takes into account the fact that the source distorts 
the medium around it. This can be achieved by making a substitution 
in In the case of a static source of mass M, this would replace its mass by 

an effective mass Q. Once this replacement is done, we find that the tree level 
amplitude is discontinuous, however, this is not problematic because simultaneously 
the region in which the expression (j^^ is applicable, i.e. r ^ r*, goes to infinity 
according to (j52j) . Hence, no vDVZ discontinuity remains in the theory. 

We also note that although the solution is nonperturbative, nevertheless the 
fields in the metric remain weak (much smaller than 1) all the way down to distances 
^ > "^M- Hence, the solution that we find never exhibits the strong field behavior 
as long as r > r^; this is similar to the conventional 4D Schwarzschild solution of 
massless gravity. 



6 Weak-field versus Cjv-expansion 

Consider small perturbations h about a flat metric t] 

Qi^u = Vfiiy + hf^u- (57) 

The weak- field expansion (WFE) is a power series expansion in /i <^ 1. In conven- 
tional Einstein's theory the WFE coincides with an expansion in Newton's coupling 
Gn- However, in the DGP model this is not so because there is another dimen- 
sionfuU parameter in the theory, mc, that contaminates the (jAr-expansion. To see 
this in some detail let us first look at the WFE and (jAr-expansion in Einstein's 
theory. This can be done in the Lagrangian or in equations of motion. We choose 
the former. Ignoring the indexes that are not important for our purposes, the WFE 
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of the Einstein Lagrangian reads: 

Ml^(hd'^h + hdhdh + +hT , (58) 

where the dots denote terms that contain only two derivatives but higher powers 
of h. The last term in (|58|) describes the interactions of gravity with matter. The 
expansion in the parenthesis of (j58p in powers of dimensionless field h is acceptable 
as long as h <^ 1 ^. One the other hand, one could rescale the field h —>■ h/M-p\. 
Then the Lagrangian takes the form 

(^d'^h + ^G^hdhdh + ..^ + ^[chihT . (59) 

The rescaled field h has the canonical dimensionality. The Newton constant emerges 
only in the graviton interaction vertices. Therefore, one can develop the standard 
Feynman diagram technique as an expansion in powers of Gat. The results of this 
expansion would coincide with the results of the WFE. 

The above arguments do not hold in general in theories where gravity gets modi- 
fied at some large distance scale Tc- This is because the new dimensionfull parameter 
Tc = enters the expansions. Below we concentrate on the DGP model to discuss 
this issue. The 4D part of the Lagrangian in the DGP model can schematically be 
written as follows: 

(hd'^h + hdhdh + ...) + m^Ml^ [hdh + hhdh + ...) + hT . (60) 

Let as look at the two terms in (j60|) . The cubic and higher powers in the first 
parenthesis can be neglected w.r.t. the quadratic terms as long as <C 1. Likewise, 
the cubic and higher powers in the second parenthesis can be neglected w.r.t. the 
quadratic terms as long a,s h <^ 1. However, the cubic term in the first parenthesis 
cannot be neglected w.r.t. the quadratic terms in the second parenthesis unless 
the derivative of the filed is very small dh -C m^. To see why this is important one 
should look at the trace equation (since we dropped the indices in all the expressions 
above we have not made a distinction between the traceless and trace Einstein 
equations) . The trace equation is subtle because the linearized bulk equations make 
the coefficient of the first quadratic term in (jHOl) vanish identically and the nonlinear 
term is the leading one Therefore, at short distances dh ^ rric some of the 

nonlinear terms in (|6U|) cannot be neglected. This is despite of the fact that the 
fields are weak (/i -C 1)! 

Let us now turn to the G Ar-expansion of (jUnjl . It is clear that the rescaling 
h —>■ h/Mpi does not lead to a Lagrangian with a single dimensionfull coupling as 
in fl59|) . Instead we get the nonlinear vertices that contain as well as rric- This 
leads to dramatic consequences. Certain nonlinear but tree-level Feynman diagrams 

■^We assume that we are in a regime of applicability of the Einstein theory itself, i.e., energy 
and momenta are smaller than Mpi 



18 



contain inverse powers of and diverge in the — limit 0. For finite the 
same diagrams give rise to the breakdown of the G 7v-expansion below the scale j2] . 
However, this is an artificial difficulty that is brought about by the expansion in Gjv- 
As we have seen in the previous sections the exact solution is regular in the irtc —>■ 
limit and fields are weak as log as r ^ vm- The metric is non-analytic in Gn showing 
that the Giv-expansion is not valid even in the regions where it naively would be 
expected to work. Nevertheless, as we discussed in the first section, the problems 
with the G]sf expansion can be fixed. Classically this is achieved by summing up the 
nonlinear corrections |2]. In the quantum theory things are more subtle since the 
summation of the diagrams is hard to perform while the expansion in G^ can lead 
to the appearance of certain higher-dimensional operators that are suppressed by 
unacceptably low scale [12]. However, as was shown in 7] and discussed in section 
1, even in this case one can formulate a quantum perturbation theory in Gn in 
which the counter-terms eliminate the dangerous loop-induced operators and the 
(jAr-expansion remains a useful tool. 

Note that a slight modification of the model can give rise to a theory in which 
the Gjv-expansion does not break down below r^, |22j. In this case the calculations 
can be performed straightforwardly as in the Einstein theory. It remains to be see 
if the approach of Ref. ,22, can be implemented in a full nonlinear theory. 

7 Existence of the bulk solution 

In the previous sections we presented an exact solution for the metric and extrinsic 
curvature on the brane. Below this will be used to argue that the solutions can be 
smoothly continued into the whole bulk. For this we use the ADM formalism |24j . 

We introduce the lapse scalar field J\f, and the shift vector field Af^ according to 
the standard rules: 

g^y = , gyy =M^ + g^.M^M'' . (61) 
After integration by parts the action (jSJ takes the form: 

S = M^ J d^xdy^-detg^^U [r + - K^^K^"") + Mj, J d^x ^f^R , (62) 

where K = g^^K^i, is the trace of the extrinsic curvature 

Kt^v = ^ (dyg^, - D^K - DM^.) , (63) 

and is a covariant derivative with the metric g^y. Note that the Gibbons- Hawking 
term implied in (jSj) is canceled in (j62j) after integration by parts. 
What we found in the previous section are the quantities: 

g^,y{x,y = 0) = g^y{x) and K^y{x,y = Q). (64) 
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This data on a brane can be considered as an initial data from which the metric 
and extrinsic curvature in the bulk could be reconstructed. Let us look at the bulk 
equations of motion that follow from by taking variations w.r.t. A/", and g^^'' 
we get respectively 

R = K'-Kl, (65) 
D^K = D'^K^,, (66) 

~g,,D^M + 2D^^{K^^M,) + g,,D\Kx) - D\NxK^,) + 2NkIK,^,^ , (67) 

where d = ^J—detg^ and G^^, denotes the Einstein tensor. We look at (jHHl -IH7 jl as 
at a system of differential equations in the y variable. Then, Eq. is just an 

algebraic equation since it contains no y derivatives except those that are already 
reabsorbed into the definition of K. The same is true for Eq. ()66j) . Furthermore, 
it is not difficult to check that (fHKjl is satisfied at y = 0+ if (j2HI)-(|SH) are fulfilled. 
Moreover, Eq. when considered as an equation defining 'yy{r,y = 0+), is 

satisfied by our solution at y = 0+. Hence, the only true differential equation which 
evolves the initial data into the bulk is ()67|1 that contains first derivative w.r.t. y on 
its l.h.s. 

Given the initial value formulation (jHK|l - (jH7jl local existence of the bulk solution 
is guaranteed |i^|2ni- The problem of global existence of the bulk solution, i.e. the 
geodesic completeness of the bulk solution obtained by the evolution equation (jU7j) . 
is not easy to establish in general. However, in the present context we can take an 
advantage of the fact that the bulk metric asymptotes to Minkowski space. Indeed, 
a number of theorems exist for asymptotically fiat spaces (see |^). In particular, in 
P7] the evolution along a timelike direction of data given on a 3 dimensional surface 
is shown to be smooth and geodesically complete under the assumption of strong 
asymptotic fiatness and a smallness condition on the initial data. The smallness 
condition is: 

sup {{dl + ifRl) + / j:idl + ir^VK^ + [ Udl + ly^'W^B' < oo , (68) 

1=0 '' 1=0 

where do is the geodesic distance from a point o in the initial data surface and B is 
the curl of the traceless part of -R^j,. Our initial data is strongly asymptotically fiat 
since 

K,, = O(r-i), 

for sufficiently large r, and satisfies the smallness condition as well. This suggests 
that the solution can be extended into a smooth, geodesically complete and asymp- 
totically fiat bulk, however, does not constitute a full proof. 
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We also mention that one could easily find a linearized form of the solution at 
r ^ Tc- This solution coincides with the linearized 5D Schwarzschild metric in which 
the integration constant is fixed to the screened mass, i.e., rj^ = tm- 

8 Concluding remarks 

In this section we discuss certain interesting issues that cinsG cis ci byproduct of our 
studies and for which further detailed work is needed. 

(i) Geodesic motion of matter and light. It is interesting to look at the geodesic 
motion of matter and light in the metric (fT2|) . The geodesic equation for the light 
rays propagating on the brane is dr/dt = e^. The useful form of A is given in Eq. 
fl2()|l . From this we conclude that for r ^ r^, the light propagates on the brane as it 
would in conventional 4D Einstein gravity with some corrections that are negligible 
everywhere except in the region r ~ r*. As we discussed, for a solar mass object this 
distance is at 10^° cm and the above effect will be overshadowed by gravitational 
effects of other sources. 

The matter geodesies are more complicated however. The 5D geodesic equation 
on the brane contains derivatives of 7 across the brane that multiply dy. Since 
the 7y is singular on the brane, it can give rise to finite contributions even though 
it is multiplied by the vanishing differentials. If so, the motion of 5D matter will 
get additional contributions from the extrinsic curvature of the brane. Also we 
could not manage to solve the matter geodesic equations, it looks like that the 
extrinsic curvature part is canceling, according to (@)), the modifications of 4D gravity 
that appear in A. On the other hand, 4D geodesies for localized matter in 4D are 
determined by the 4D induced metric governed by A. 

(a) Comments on Black Holes. In the present work we were dealing with a 
macroscopic source on the brane, a star for instance. In the context of the DGP 
model these type of sources are assumed to be localized on a brane by a certain 
mechanism not related to gravity itself. If the sources were not localized, then, as is 
known [2H|, the brane with the induced graviton kinetic term has effectively repulsive 
gravity and it would push any source off the brane. For instance, ordinary black 
holes cannot be held on the brane. However, charged black holes could still be quasi- 
localized if the corresponding gauge fields are localized. It would be interesting to 
see how the properties of charged black holes would differ. It is also very interesting 
to understand in detail the structure of the horizon of a black hole in the bulk. 
Our preliminary findings suggest that it should have a cylindrical form, where the 
cylinder extends from a brane into the bulk to a distance that is bigger that tm 
but smaller than r^:. Further detailed investigations are needed to understand the 
validity and implications of these observations. 
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